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ABSTRACT 

Context. Cosmic rays are thought to be accelerated at supernova remnant (SNR) shocks, but conclusive evidence is lacking. 
Aims. New data from ground-based y-ray telescopes and the Large Area Telescope on the Fermi Gamma-ray Space Telescope are 
used to test this hypothesis. A simple model for y-ray production efficiency is compared with measured y-ray luminosities of SNRs, 
and the GeV to TeV fluxes ratios of SNRs are examined for correlations with SNR ages. 

Methods. The supernova explosion is modeled as an expanding spherical shell of material that sweeps up matter from the surrounding 
interstellar medium (ISM). The accumulated kinetic energy of the shell, which provides the energy available for nonthermal particle 
acceleration, changes when matter is swept up from the ISM and the SNR shell decelerates. A fraction of this energy is assumed to be 
converted into the energy of cosmic-ray electrons or protons. Three different particle radiation processes — nuclear pion-production 
interactions, nonthermal electron bremsstrahlung, and Compton scattering — are considered. 

Results. The efficiencies for y-ray production from SNRs for these three processes were compared with y-ray luminosities from 
Tycho, RX J 17 1 3. 7-3946, and other SNRs. SNRs may suddenly lose y-ray power at > 10 4 yr. Our results are consistent with the 
hypothesis that supernova remnants accelerate cosmic rays with an efficiency of « 10% for the dissipation of kinetic energy into 
nonthermal cosmic rays. Weak evidence for an increasing GeV to TeV flux ratio with SNR age is found. 
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1. Introduction 

Supernovae have been hypothesized (Ginzburg & Syravotskii 
1964 Hayakawa 1969) to be the sources of the cosmic rays, 
which consist primarily of relativistic protons, ions and elec- 
trons. Firm evidence for this hypothesis is lacking because cos- 
mic rays are deflected by intervening magnetic fields and do 
not trace back to their sources. A distinctive y-ray signature of 
cosmic rays results from interactions of cosmic-ray protons and 
ions with gas and dust in the vicinity of their acceleration sites, 
making y rays with a characteristic n° — > 2y emission feature 
peaking at 67.5 MeV in a photon number representation (Stecker 
fT97TT >. 

Resolution and sensitivity limitations of previous y-ray tele- 
scopes made it impossible to precisely localize y-ray sources in 
the Galaxy and to measure the expected 7r°-decay spectral sig- 
nature. This has changed with the new generation of y-ray tele- 
scopes, including AGILE (Astro-rivelatore Gamma a Immagini 
LEggero; Gamma-ray Light Imaging Detector) and the Fermi 
Large Area Telescope (LAT). Preliminary evidence for a cut- 
off below several hundreds of MeV from some SNRs has been 
presented (Giuliani et al. [20TTl Abdo et al. 1201 Pel l2TjTTjel ) as evi- 
dence for this feature, so this one-hundred year old problem may 
soon be solved. 

While spectral analyses at photon energies E y < 100 MeV 
are still in progress, we can also ask whether the measured y-ray 
luminosities of SNRs are consistent with expectations of effi- 
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ciency for cosmic -ray acceleration at SNR shocks. Furthermore, 
trends in the y-ray spectra might reveal a progressive behavior 
characteristic of cosmic-ray acceleration in these sources. 

Here we test the hypothesis that cosmic rays are accelerated 
at SNRs by calculating the fraction of supernova kinetic energy 
that must be converted into cosmic-ray particles in order to pro- 
duce the y rays that are observed from some of these remnants. 
We also plot the GeV to TeV energy-flux ratio with age. 

2. Particle acceleration at SNR shocks 

Supernovae (SNe) are found in two main classes, namely ther- 
monuclear Type la SNe and core-collapse Type II SNeQ The 
progenitors of Type la SNe are Solar mass white dwarfs com- 
posed of carbon and oxygen. Thermonuclear burning of the ma- 
terial in the white dwarf releases a prodigious amount of energy 
and creates a shock wave as the burned material expands. For 
stars more than w 8 times as massive as the Sun, in contrast, the 
stellar interior burns until much of it has been converted to an 
Fe core, which will collapse when electron-degeneracy pressure 
is no longer sufficient to support the core against gravity. A core 
bounce sends a shock wave out that expels a portion of the stel- 
lar envelope (for a general review, see Burrows 2000). Neutron 
stars can be formed in Type II SNe, but not in Type la SNe. 

In both Types la and II SNe, the energy available for nonther- 
mal particle production is derived from the directed bulk kinetic 

1 Other subclasses of core-collapse SNe subclasses are Type lb and 
Ic SNe. 
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energy of the shell. This energy is dissipated in the form of the 
random kinetic energy of matter swept up at the SNR shock. 
The amount of energy that is extracted and transformed into 
cosmic-ray energy is highly uncertan, and represents a major 
open question in SNR modeling. Different prescriptions related 
to the number of swept-up particles or the kinetic energy dissi- 
pated at the shock front have been considered (e.g., Sturner et al. 
[19971 Reynolds [19981 Baring et al. [19991 Tang et al. l20TTl) . but 
the exact conversion efficiency depends on poorly understood 
microphysical processes (e.g., Blasi et al. [20051 . Here we re- 
consider this problem using an idealized model of an expanding 
spherical shell of matter and a uniform surrounding medium (cf. 
Truelove & McKee lT999l Finke & Dermer l20T2l l. 

Let Eo = 10 51 £ , 5i erg represent the total kinetic energy of the 
supernova ejecta with mass Mo = mM Q , where M Q = 2.0 x 10 33 
g is the mass of the Sun. For Type la and II SNe, m = 1 .4 and 
m = 10, respectively, with £51 =s 1 for both types. The SN is 
here assumed to explode in a medium of uniform density no 
and form an outward moving spherical shell of material with 
radius r = r(t) at time t after the explosion. As it expands, the 
shell sweeps up material with mass M su (t) = Anr 3 po/3, where 
the density po = tn p na, and the ISM is assumed to be com- 
posed of hydrogen. The initial speed of the shell is therefore 
v = V2£o/Afn = 10 4 V^si/mkms" 1 . 

The Sedov radius rs e d is defined by the condition that the 
swept-up mass equals the explosion mass, that is, Anp$rl J3 = 

M , implying r Sed = (3M /4np ) 1/3 = 6.6 x I0 n (m/n ) h3 cm. 
The Sedov time tsed is the time when the shell transitions from 
the coasting phase at r < rs e d to the Sedov phase at r > r^ed, and 
is given by t Sed = r Sed /v = 210 m 5/6 E~l /2 n~ 1/3 yr. 

For nonrelativistic speeds v„ <K c and adiabatic shock waves, 
the kinetic energy 



E ke = \[M Q + M su (t)}v 2 



(1) 



is conserved, so that dEk e {t)ldt = 0. Introducing dimensionless 
radius x = r/rs e d and time t = f/fsed in units of Sedov radius 
and Sedov time, respectively, and noting that v = dr/dt, Eq. (Q]i 
leads to the equation of motion 



(1 + x 3 ) l/2 dx = dr 



(2) 



for the dimensionless shell radius x. From this, the asymptotes 
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Fig. 1. Top: Red curve labeled "1" shows the rate at which ki- 
netic energy is swept into the SNR shell, and the blue curve la- 
beled "2" shows the rate at which kinetic energy is lost due to de- 
celeration of the shell. The net rate of change of kinetic energy in 
the SNR shell is given by the green curve, labeled "3". Bottom: 
Integrated kinetic energy corresponding to the three terms shown 
in the upper panel. Curve 3 gives the accumulated swept-up ki- 
netic energy, E su (f). 
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and 

V(T) 



Vo , :«1 
v (5r/2)- 3/5 , t » 1 



(4) 



can be derived. 

If nonthermal particles extract their energy from the kinetic 
energy of swept-up matter, then the change of kinetic energy of 
swept-up matter with time is given by 



dE™(t) _ d_ 
dt dt 



1 



:M.„ 



1 



+ M,„vv 



„ 21 An , 

2nr v p + — r vvp . 



(5) 



The first term on the right-hand-side of Eq. (|5]l is the rate at 
which kinetic energy is swept into the shocked material, and the 



second term is the rate at which the kinetic energy of the swept- 
up matter is lost by deceleration of the shell (and is therefore 
negative). The sum represents the net rate of change of kinetic 
energy in the SNR shell. These terms are are shown in the top 
panel of Fig. [T] and the integrated energy associated with each 
term is shown in the bottom panel of Fig. Q] Here we let £51 = 1, 
m = 1.4, and ran = 0.1 cirT 3 . Integrating the swept-up kinetic 
energy leads to a logarithmically diverging total energy that ex- 
ceeds the initial kinetic energy of the explosion. By adding the 
change of energy associated with deceleration, the summed con- 
tribution recovers the initial kinetic energy of the explosion at 
late times. The process of sweeping up matter from the surround- 
ing medium has transformed the directed kinetic energy into in- 
ternal thermal and nonthermal kinetic energy of swept-up matter. 



Note how rapidly, > 



-2.5 



, that the total injection rate of kinetic 



energy decreases in the Sedov phase. 
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3. Cosmic-ray interactions in the SNR environment 

Cosmic-ray ions and electrons can interact with target parti- 
cles to make y rays through secondary nuclear production and 
bremsstrahlung, respectively. Additionally, relativistic electrons 
will Compton-scatter photons of the ambient radiation fields to 
the y-ray regime (see Blumenthal & Gould 119701 and Dermer 
& Menon 2009 for reviews of high-energy astrophysical radi- 
ation processes). There are three targets available for particle- 
particle interactions, namely from (1) the explosion mass, (2), 
the swept-up mass, and (3) the ISM gas. The explosion mass 
density n ex {t) = Mo/m p V s h(t), where the volume of the expand- 
ing shell of width A is V sh (t) = 47r[r 3 - (r - A) 3 ]/3 -> 47rr 2 A in 
the limit u = A/r «: 1 . Hence 



«ex(0 



no 
3x 3 u 



It, 

3u 1 (5t/2)- 6 / 5 



T«l 

r » 1 



(6) 



The density of the swept-up mass is increased by the compres- 
sion ratio. For a strong shock, the shell of swept-up matter has 
density n m = 4«n. The most efficient radiation production there- 
fore occurs when the cosmic rays interact with explosion mass 
in the coasting phase, and swept-up mass in the Sedov phase. 
These two regimes can be bridged with a formula for the maxi- 
mum target density available for interactions, given by 



^mflA'(0 — 



no 

3 MX-' 



(1 + 12ux 3 ) . 



(7) 



Secondary nuclear production makes y rays from the decay 
of neutral pions, as well as through emissions from pion-decay 
electrons and positrons. For cosmic-ray protons with energies 
E p = m p c 2 y p well above the pion-production threshold energy 
at E p = 300 MeV, the rate at which the proton Lorentz factor 
y p changes is given by -y pp = K p ccr pp n(t)y p , where K p w 0.5 
is the mean inelasticity (fraction of energy lost per collision), 
and o-pp x 30 mb is the p-p inelastic cross section, so that the 
energy-loss timescale associated with secondary production is 



'pp 



\lplj. 



pp\ 



[K p ccrppn(t)]- 1 s 2.2 x 10 ls /n(0 s. The max- 
imum y-ray luminosity that can be made by secondary nuclear 
production is therefore 



L„p(t) = 



rip E su {t) 
lt P p(t) ' 



(8) 



with n(t) replaced by n max (t), from Eq. in t pp (t). The factor 
1/3 represents the mean fraction of energy lost in an inelastic 
nuclear interaction that emerges as y rays, and the term t] p is the 
fraction of accumulated swept-up kinetic energy that is trans- 
formed into cosmic-ray protons. 

The bremsstrahlung energy-loss rate for electrons with en- 
ergy E e = m e c 2 y e interacting with ions of density nz and 
charge Z in a fully ionized medium is given by -y// = 
kffafCO-j[EnzZ(Z+l)]y t ,, where a/ and cr T are the fine-structure 
constant and Thomson cross section, respectively. The term 
k ff = (3/2n)(ln2y e - |) * 4.3, taking y e * 6000 for 
the typical electron Lorentz factor y e that makes GeV y rays 
through this process. The characteristic bremsstrahlung energy- 
loss timescale is therefore ?//(?) = \y e ljff \ — 1-6 x 10 15 /n(f) s, 
and the maximum bremsstrahlung luminosity is 



L f f{t) = 



TleJfEsuif) 
tffit) ' 



(9) 



again replacing n(t) with n max (t) in f//(f)- 

The loss rate of cosmic-ray electrons by Compton scattering 
and synchrotron radiation is given by -yes - (4/3)c<x x (t/ r + 



UB)y 2 /m e c 2 , where U y and Ub are the energy densities of 
the photon and magnetic fields, respectively. The characteristic 
energy-loss timescale through Compton and synchrotron pro- 
cesses is therefore ?cs(s)= \y e /jcs\ = 7.7 x 10 19 /y<,(l + q) 
« 7.1 x 10 13 /(1 + q) V^GeV, where q = Us/Uy, and we assume 
that the dominant radiation field is the CMBR, with energy den- 
sity Uy = Ucmb — 4.0 x 10 13 erg cm -3 . In the final expression, 
y e is replaced by the value of the electron Lorentz factor that will 
scatter a typical CMBR photon to a y ray with energy Ecev GeV. 
The y-ray luminosity from Compton scattering is therefore 



L c (t) = 



t] e ,c E su (t) 

qtcs(t) 



(10) 



because only the fraction q~ l of the total nonthermal energy is 
radiated as Compton-scattered y rays. 

The curves shown in Fig. [2] give results for the energy-loss 
timescales in the top panel, and the y-ray luminosities of SNRs 
in the bottom panelQ For the model results, we assume that Type 
la and II SNe involve 1.4 and 10 M of explosion mass, and that 
the surrounding medium densities are 0.1 and 10 crrT 3 , respec- 
tively. The total kinetic energy of the explosion is assumed to be 
10 51 erg in both cases, and we let u = 0. 1 for the relative width of 
the shell to the radius. The Sedov timescales for these parameters 
are 595 y and 660 y for the Type la and Type II SNRs, respec- 
tively. Also shown are the energy-loss timescales for cosmic-ray 
electrons that scatter CMB photons to y rays of energy E 7 . Here 
a magnetic field energy density 10 times as large as the photon 
energy density (q = 10) is assumed. The solid line in the up- 
per panel separates the weak and strong cooling regimes; in the 
strong-cooling regime our treatment does not apply. 

The lower panel of Fig.|2]gives model results for the y-ray lu- 
minosities corresponding to the same parameters as shown in the 
upper panel. The swept-up kinetic energy is assumed to be trans- 
formed into cosmic-ray protons and electrons with efficiencies 
of rjpp = rjejf = 10%, and r\ e ,cs = 0.1% for the three processes. 
The low efficiency for converting kinetic energy into electrons 
is small for the latter process, but these electrons have to have 
very large Lorentz factors with y e = 10 6 V^GeV- Note that the 
increase in swept-up kinetic energy competes with the decrease 
in target density so that the y-ray luminosity is roughly constant 
with time for secondary nuclear production and bremsstrahlung. 
By comparison, the target photon density for Compton scattering 
is constant, so the y-ray luminosity through this process tracks 
the increasing swept-up kinetic energy shown in Fig.Q] 



4. Comparison with data 

We calculated the > 300 MeV y-ray luminosities L y from the 
data provided in the references for the SNRs listed in Table 1 . 
The large error bars in L y reflect both the uncertainties in the 
y-ray measurements as well as distance uncertainties. The val- 
ues of L y range from w 10 34 erg s~' to » 3 x 10 36 erg s~'. 
The data show an apparent increase in luminosity from young 
to middle-aged SNRs, followed by a decline for the older SNRs, 
but a larger sample of L y from SNRs will be needed to confirm 
this trend. Most of these SNRs are Type II, so that if the pa- 
rameters used to characterize the Type II SNRs are reasonable, 
then we find that 77 ~ 10% for nuclear and bremsstrahlung pro- 
cesses, consistent with expectations of production efficiency if 
SNRs accelerate the cosmic rays. Also note that Type la SNe are 



2 We calculate a y-ray luminosity for CTB 37A larger, by an order of 
magnitude, than Castro and Slane (2010ll. 
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Fig. 2. Top: Energy-loss timescales of cosmic -ray protons 
through secondary nuclear production (t pp ; solid curves) and 
of cosmic-ray electrons through bremsstrahlung (tff, dashed 
curves) are shown for parameters corresponding to Type la 
and Type II SNe, as labeled. Also shown are time scales for 
Compton and synchrotron energy losses (fc; dotted curves) of 
cosmic -ray electrons that scatter CMB photons to y rays of en- 
ergy E y . Bottom: Theoretical estimates of y-ray luminosities of 
SNRs through secondary nuclear production (L pp ; solid curves), 
bremsstrahlung (Lff, dashed curves), and Compton scattering 
(Lc; dotted curves) are shown by curves. Data points show the 
measured > 300 MeV y-ray luminosities for the SNRs listed in 
Table 1. 



predicted to be m 2 orders of magnitude less y-ray luminous than 
Type II SNe, so that there should be a preference to detect this 
latter class, in addition to the greater rate of Type II SNe in the 
Galaxy. Very small efficiencies are needed to make y rays from 
Compton scattering, but except in the case of RX J17 13.7-3946, 
the broadband y-ray spectra of these SNRs disfavors a Compton- 
scattering origin. 

We also plotted the ratio of the 0.3-3 GeV energy flux to the 
300 GeV - 3 TeV energy flux for the SNRs listed in Table 1 . The 
result is shown in Fig.|f] All of the SNRs except for RX J 17 13.7- 
3946 have GeV to TeV flux ratios of order unity or greater. The 
middle-aged SNRs are more luminous in GeV y rays than in TeV 
y rays by factors of > 10 - 100. RX J17 13.7-3946 is anoma- 
lously low in comparison. There is a weak trend for increasing 
GeV to TeV flux ratio with age, though if RX J1713.7-3946 is 
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Table 1. GeV-TeV SNRs. 



SNR 


Type Age (kyr) d (kpc) 


Refs. a 


Cas A 


II 


0.33 


3 4 +03 


1 


Tycho 


la 


0.44 


1.5-4.0 


2 


RX J1713.7 


-3946 H(?) 


1.6 


0.9- 1.7 


3,4 


CTB 37A 


II 


2 ±0.5 


10.3 + 3.5 


5,6 


W49B 


II 


1-4 


8-11 


7 


W51C 


11 


10-30 


7 ±1.5 


8 


W44 


11 


15-20 


~ 3 


9, 10 


W28 


II 


35-45 


1.9 ±0.3 


11, 12 


IC443 


Ib(?) 


20 - 30 b 


~ 1.5 


13 



a 1. Abdo et al. l2010al 2. Giorda no et a l. l20T2l 3. Abdo et al.l2TJTTl 

4. Fesen et al. [20121 5. Brandt |20TT] 6. Aharonian et al.l2008l 7. 

Abdo et al. l2010bl 8 Abdo et al. |2009l 9. A bdo et al. l2010cl 10. 

Wolszczan et al. [19911 11. A bdo et al. l2010dl 12. V elazquez et al. 

[20021 13. Abdo e t al. l2010el see also Caprioli|20TT] 
b Troja et al. 2008 argue from XMM-Newton X-ray observations that 

the age is < 10 kyr. 
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Fig. 3. Ratio of 0.3 - 3 GeV to 0.3 - 3 TeV fluxes for the SNRs 
given in Table 1 . 



removed, the trend is tentative at best. This might be expected 
if the dominant processes are Compton scattering and nuclear 
production, noting from Fig.[2]that cosmic -ray electrons making 
TeV photons lose their energies on timescales of ~ 10 4 y, leav- 
ing the long-lived cosmic -ray protons that would preferentially 
make GeV radiation. 

Compton scattering as a TeV y-ray emission process for 
middle-aged SNRs might be ruled out by noting how rapidly 
the high-energy electrons that scatter CMBR photons to TeV en- 
ergies cool (Fig. [2]>. If these electrons are powered by swept-up 
kinetic energy, the rapid decrease in the kinetic-energy injection 
rate might be inadequate to replenish these electrons, as shown 
by curve 3 in Fig. [TJ Indeed, a rapid decline in the X-ray syn- 
chrotron radiation from y ~ 10 8 / ^J(BJToJtG) electrons shortly 
after the Sedov phase is found by plotting the nonthermal X- 
ray luminosity as a function of radius (Nakamura et al. 120121 1. 
which is closely related to the remnant age. Reduction in the 
GeV - TeV y-ray luminosity in middle-aged SNRs, if made by 
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primarily by protons and ions, seems to mirror the behavior of 
X-ray synchrotron emission from high-energy electrons at ear- 
lier times. The decline in the y-ray luminosity would, however, 
instead be due to escape rather than energy loss, given the long 
timescale for energy loss of cosmic -ray protons through inelastic 
nuclear processes. Finally, we note that the assumption of a uni- 
form medium surrounding the SNR is an extreme simplification 
of the actual environment of SNRs, particularly the core-collapse 
type found in star-forming regions. 
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5. Conclusions 

Data from the new generation of y-ray telescopes can be used 
to test the hypothesis that SNRs are the sources of the cosmic 
rays. We have constructed a simple model for y-ray production 
in SNRs, and find that for parameters that are expected to ap- 
ply to Type la and II SNe, « 10% conversion efficiencies of di- 
rected kinetic energy to nonthermal protons or electrons are re- 
quired to account for the measured y-ray luminosities of SNRs 
if the radiation mechanism is secondary nuclear production or 
bremsstrahlung, respectively. Smaller efficiencies can account 
for the y rays if the emission process is Compton scattering, but 
the energy has to be converted to highly relativstic electrons. We 
also looked for trends in L y and the GeV to TeV flux ratios with 
SNR age. SNRs are found to be most y-ray luminous at ages 
of a: 10 3 - 10 4 y. The decline of L y at later times might arise 
from escape of cosmic rays from the acceleration and dense tar- 
get interaction sites. The GeV to TeV flux ratio of SNRs also 
displays a tendency to increase with age, but additional analyses 
and study are needed to establish these trends. 
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